Variational procedure for nuclear shell-model calculations and energy-variance 
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We discuss a variational calculation for nuclear shell-model calculations and propose a new pro- 
cedure for the energy-variance extrapolation (EVE) method using a sequence of the approximated 
wave functions obtained by the variational calculation. The wave functions are described as linear 
combinations of the parity, angular-momentum projected Slater determinants, the energy of which is 
minimized by the conjugate gradient method obeying the variational principle. The EVE generally 
works well using the wave functions, but we found some difficult cases where the EVE gives a poor 
estimation. We discuss the origin of the poor estimation concerning shape coexistence. We found 
that the appropriate reordering of the Slater determinants allows us to overcome this difficulty and 
to reduce the uncertainty of the extrapolation. 

PACS numbers: 21.60.Cs, 27.50.-|-e, 24.10. Cn 
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I. INTRODUCTION 

The role of large-scale shell-model calculation has been 
increasing in nuclear structure physics with the recent de- 
velopment of faster parallel-computing capabilities. The 
most popular method to perform shell-model calculations 
is the Lanczos method and its variants, which can han- 
dle 0(10^^) configurations. However, the feasibility of 
this method is still hampered by the exponential increase 
of the Hilbcrt space of the model space as a function 
of the number of nucleons. To overcome this difficulty, 
much effort has been given to developing approximation 
schemes to an exact shell-model diagonalization method 
[J-Q. These approximation methods truncate the whole 
Hilbert space to a relatively small subspace which is de- 
termined by various sophisticated methods. Therefore, 
the approximated energy in the subspace always provides 
us with an upper limit to the exact energy because of the 
variational principle, and an unavoidable small gap re- 
mains between the approximated energy and the exact 
energy. 

This gap can be removed by extrapolation, which has 
been intensively studied from various perspectives 
The fundamental ingredient of these studies is the ex- 
trapolation of the approximated eigenenergy by expand- 
ing the subspace into full Hilbert space. One realiza- 
tion of this spirit is the Exponential Convergence Method 
(ECM) d, in which the approximated eigenenergy is ex- 
trapolated as an exponential function of the dimension 
of the truncated subspace. Another scheme of the ex- 
trapolation is utilizing the energy variance of a sequence 
of the approximated wave functions fl2!| . In this scheme. 



the approximated eigenenergy is extrapolated as a poly- 
nomial function of the corresponding energy variance be- 
cause the energy variance of the exact wave function van- 
ishes. The uncertainty of the extrapolated energy is ex- 
pected to be small because the function used for extrap- 
olation is a first- or second-order polynomial, while the 
corresponding function is exponential in the ECM. This 
energy-variance extrapolation (EVE) was introduced in 
condensed matter physics [l8| and applied to various ap - 
proximations in nuclear shell-model calculations [l^, Ej- 



In this article, we report how the EVE method works 
with the approximated wave function represented by the 
linear combination of the parity, angular- momentum pro- 
jected Slater determinants, such as the Monte Carlo Shell 
Model (MCSM) [B]- A case in which the EVE is diffi- 
cult is reported in Ref. We discuss the origin of 
this case with relation to shape coexistence and how the 
case is solved by choosing an appropriate reordering of 
the Slater determinants. 



We explain the variational calculation to generate a se- 
quence of approximated wave functions using the MCSM 
and conjugate gradient (CG) method ^gOj in Sect|lTl In 
Sect lllli we discuss how the simple EVE method works. 
We also show a typical case in which the EVE shows 
a poor result. In Sect II VI we introduce the reordering 
technique of the extrapolation method to remedy this 
difficulty, and demonstrate the feasibility of the method. 
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II. APPROXIMATED WAVE FUNCTIONS 

We briefly describe how to construct the truncated 
subspace to the full Hilbert space in the nuclear shell- 
model calculations. The approximated wave function is 
written as a linear combination of angular-momentum- 
projected, parity-projected Slater determinants, 

N J 

l*^)-E E fnJ<PirK\^n). (1) 
n=l K=-J 

where A'^ is the number of the basis states, which corre- 
sponds to the "MCSM dimension" in Ref.[2l|. is 
the angular-momentum, parity projector defined as 

(2) 

where = (a, /3, 7) represents the Euler angles, and 
P^MK denotes Wigncr's I?-function. 11 denotes the parity 
transformation. Each |0„) is a deformed Slater determi- 
nant, 

\^n)^X{{^D\t>c\^\^), (3) 

where c| denotes a creation operator of single particle 
state i and |— ) denotes an inert core. The coefficient f^x 
is determined by the diagonalization of the Hamiltonian 
matrix in the subspace spanned by the projected Slater 
determinants, -PMi^l*^")- This diagonalization also deter- 
mines the energy, Em = n\H\^ n) , as a function of TV. 
Note that the dimension of the subspace is the product 
of the number of states, iV, and the degree of freedom of 
the z-component of angular momentum, 2 J -|- 1. We in- 
crease N until Em converges enough, or the extrapolated 
energy converges. 

The coefficient D*^") is given by a variational calcula- 
tion utilizing the auxiliary field Monte Carlo technique 
and the CG method. Each D*^"^ is determined by min- 
imizing 

EN=n while keeping D^^)^ _fj(n-i)^ 

m a 

manner similar to the MCSM [ij, the few-dimensional 
basis approximation and the Hybrid Multidetermi- 
nant method @. We perform the MCSM procedure in 
a small number of steps to get the initial states of the 
CG process in order to avoid a trap by local minima. 
Since the D^") is determined sequentially, hereafter, we 
call this procedure the Sequential Conjugate Gradient 
(SCG) method. 

The energy variance of the approximated wave func- 
tion \^m) is also evaluated as {AH'^)n = {^n\H'^\'^n)- 
Eff. The energy, energy gradient, and energy variance 
are evaluated under the angular-momentum, parity pro- 
jection technique throughout this work. In this sense, 
our method is a "variation after projection" . 



III. ANOMALOUS KINK IN THE EV PLOT 

The ^^Ge is a typical case such that the relation be- 
tween the energy and its variance is not monotonic, which 
is reported in Ref. [l^ . This nucleus exhibits a feature of 
shape coexistence due to the A'' = 40 magicity [2^, and 
this is considered to be a main reason for ill-behavior as 
will be discussed later. In the ^^Ge shell-model calcu- 
lation, we take the /sp^g-shell, which consists of O/5/2, 
1P3/2J ^Pi/2i and 0(7g/2 single-particle orbits, as a model 
space and the effective interaction JUN45 is used [2^ . 
The TO-scheme dimension of ^^Ge reaches 140,050,484, 
which can be handled directly by recent shell-model di- 
agonalization codes. The exact values in Fig. [ijrepresent 
the results of the shell-model diagonalization utilizing the 
code MSHELL64 

Figure [1] shows the energy and energy variance of j^'^v) 
with 1 < TV < 100 for the Of state. The I'S/n) is obtained 
by the SCG method. Hereafter, we call the plot of the 
energy as a function of a variance "EV plot". In the 
EV plot concerning the 0^ state, one can find a similar 
anomalous kink in the case of the MCSM [19]. This kink 
reduces the region available for the second order fit, and 
deteriorates the certainty of the extrapolation, resulting 
in a 50 keV overestimation. To obtain the Oj energy, 
additional 100 bases are generated by the SCG method 
for minimizing iS^vCOj). The dashed line in Fig. [1] is 
fitted for the points concerning J = Oj, whose extrap- 
olation agrees with the exact energy well, but a small 
underestimation remains. 




Energy Variance (MeV ) 

FIG. 1: (Color online) Energy vs. energy variance plot of the 
ground state and state of '^^Ge by the SCG method. The 
circles and triangles show the energies and energy variances of 
0^ and O2, respectively. The lines are fitted for these points 
using a second-order polynomial. 

Figure [2] shows the total energy surface provided by 
the Q-constrained Hartree-Fock calculation [l^] using 
the same shell-model Hamiltonian. There are two low- 
energy regions corresponding to the shape coexistence 
phenomenon [2^ . In order to discuss the intrinsic struc- 
tures of the of and Oj wave functions, we plot the de- 
formation of each unprojected basis state of the 
SCG wave function I^'at) in Fig. [5J The location of 
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the scattered circles shows the quadrupole deformation, 
namely, {4'n\Qo\4>n) and ((/)„|Q2|</>n) where Qm is an M- 
component of the mass quadrupole operator and is 
rotated so that {(t>n\Q±i\<t>n) — 0. It is interesting to see 
that the circles scatter in a broad region of the energy sur- 
face, not near the local minima, but on the hillside and 
triaxially deformed regions due to the effect of the con- 
figuration mixing and the "variation after projection" . 
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FIG. 2: (Color online) Total energy surface of the ^^Ge by the 
Q-constrained Hartree-Fock calculation. The scattered circles 
denote the deformation of the basis states of the SCG wave 
functions for (a) 0^ ground state, and (b) second 0^ state. 
The contour width is 250 keV. 

The area of each scattered circle on Figs. [2] (a) and (b) 
is proportional to the overlap probability of each pro- 
jected basis state with the resulting many basis states, 
n\P'' \'t'n)\^ , where TV is the normalization factor, 
TV — {(j)n\P''' \4'n) ■ Concerning the 0^ state shown in Fig. 
[2] (a), the overlap probability is rather small and 0.54 
at most. The circles with relatively large overlap scat- 
ters in a broad region of the energy surface. The overlap 
probability with P'^j^i), e.g. the Hartree-Fock solution 
with the variation after projection is 0.32, which is mod- 
est. This overlap implies that the 0^ state is described 
by a linear combination of a relatively large number of 



basis states essentially, or the effect of the configuration 
mixing plays an important role. On the other hand, in 
Fig. [2] (b), the points of the large overlap concentrate 
near the spherical region concerning the 0^ state. The 
overlap probabilities corresponding to these points are 
large with the Oj many-basis state, and the probability 
is 0.67 at most. This means that the wave function can be 
well approximated by a few number of projected Slater 
determinants, or mean-field description works far better 
than in the case of the Of state. This property makes 
the first few SCG basis states for minimizing Of energy 
dominated by the Of state, not by the Of state. These 
different properties of the Of and Of wave functions give 
rise to the different behavior of energy convergence as 
a function of the number of basis states, which makes 
extrapolation with a number of basis states difficult . 

In Fig. [U the variance shows the local minimum at 
E — —182 MeV, which is near the energy of the Of state. 
This figure indicates that the wave function comprised of 
the first 10 bases are dominated by the true Oj state, not 
the of state. This is consistent with the discussion using 
the energy surface in the previous paragraph. 

Such an anomalous kink in the EV plot is also seen in 
the 8f state of ^^Ni with the FPD6 effective interaction 
[25} . ^^Ni is known to have shape coexistence [2^ , which 
is consistent with the previous discussion. This situa- 
tion may occur in the case where the next lowest energy 
eigenvalue is close to the target one and the mean-field 
solution favors the next lowest state. 



IV. EVE AND REORDERING OF BASIS 
STATES 

The emergence of the anomalous kink discussed in 
Sect lllll gives rise to a large uncertainty of the EVE. In 
order to remove such a kink and to improve the preci- 
sion, we introduce the reordering of basis states of the 
SCG wave function, represented in Eq.dl]). 

The SCG wave function comprises a set of Nm basis 
states, |<^ri) with 1 < n < Nm- These basis states also 
give us a sequence of approximated wave functions, Y^n) 
with 1 < < Nm, which is used for the simple EVE 
method. Meanwhile, reordering these Nm basis states 
using a permutation, cr(7i), yields another sequence of 
approximated wave functions. 



N 

n=l K=~J 



fnPPim\^l:°^) (4) 



where \(l>n°^) = |0cr(n))- The corresponding energy 
and the energy variance (Ai/^)^"-* provide us with a new 
EV plot. By reordering the basis states with o'(n), we can 
simulate another truncation scheme. Because the behav- 
ior of the fitted line in the EV plot depends on the trun- 
cation scheme [T3| . we remove an anomalous behavior of 
the EV plot and make the extrapolation stable by using 
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an appropriate a{n). We describe the way how to obtain 
an appropriate order of basis states in this section. 



A. Procedure of reordering technique 

The relation of the energy and its variance is usually 
assumed to be expressed as a second-order polynomial. 
Because a second-order-term error is roughly estimated 
as <5c2((A_ff^) at)^, where 6c2 is a second-order-term error 
of x-square fitting, the uncertainty of the second order 
term often causes a relatively large error of the extrap- 
olated value in case (A7?^)jv, is not small enough. If 
we can find an order of the basis states in which the 
fitted curve is a first-order polynomial, the error of the 
extrapolation will come mainly from the coefhcient of 
the first-order term, which can be estimated as roughly 
proportional to dci ( AiJ ) , where Sci is the first-order- 
term error of x-square fitting. This error should be far 
smaller than that of second-order polynomial if {AH'^)n 
is large. Thus, our strategy is to select the order of the 
basis states, in which the fitted curve is close to linear, 
by changing the order of the basis states. 

Now, we discuss the relation between energy difference 
and energy variance following the idea of Ref.[l^. We 
define the energy difference SE between the energy ex- 
pectation value E — {'i/\H\'^) of the approximated wave 
function l^t) and the lowest exact energy eigenvalue Eq 
as 



By eliminating (P, we obtain 



SE ^ mH\-^)-Eo, 



(5) 



and the energy variance of the approximated wave func- 
tion as 



(6) 



An approximate ground state \^!) can be decomposed 
into the exact eigenstate, \tpo), and the rest of the com- 
ponent, \ipr) as 



I*) = c\^Jq) +d\ijr), 



(7) 



where + = 1. iV'r) is expanded by the exact excited 
states such as 



IV'r) = ^ C„|-0n). 

By defining the moments Dj as 

D,=Y,cl{En~Eo) 



we obtain 



6E = d^Di, 



(8) 



(9) 



(10) 



(11) 



(AH^) 



SE - {SEf 



(12) 



In Eq. PT|) . {AH"^) is written as a second-order polyno- 
mial of SE, which explains the parabola shape of the EV 
plot in Fig. [T] (circles below -182MeV). The authors of 
Ref.Ts^ solve (5^; as a function of (AiJ^), which is shown 
to be approximated as a second-order polynomial. Now, 
we assume SE is proportional to (AiJ^), which is real- 



ized in case -S^- 



oo. Because of 



d(SE) 



D2 ' 

we reorder the basis states so that q^^^i ^ is as small as 
possible, resulting in linear proportionality. 

In practice, we perform the following procedure: 

1. A fixed number, N„i, of the basis states is obtained 
by the SCG method. 

2. Choose the A^,„-th basis state 1'/'^°'') from Nm can- 



didates by minimizing 



^(ro) _ (ro) 



3. Set N as — 1. Choose the iV-th basis state, 
l^^""*) from TV candidates (basis states except for 
already fixed states, with N + 1 < N' < N,n 



by minimizing 



(A//2)(^=)^_(A//2)(^°) 



4. Set iV as — 1. The previous step is iterated re- 
cursively up to determining the first state, 

Note that this procedure needs no additional heavy 
computation to the SCG, because the matrix elements of 
energy and its variances are evaluated once and stored, 
and we only require the diagonalization of the matrix 
whose dimension is A^ < A^,,, ~ 100 for each candidate 
order of the basis states. 

As a result, we obtain a fitted line which is closest to 
linear in the provided set of basis states. Because the 
reordering makes the gradient of the EV plot as small as 
possible, the anomalous kink discussed in Sect. IIIII van- 
ishes, which is demonstrated in the following subsection. 



B. 



Ge in f^pgg-sheU 



We apply the reordering technique to the ^^Ge in 
f^pgQ-sheW. The filled circles in Fig. [3] shows the EV 
plot with the reordering technique. The anomalous kink 
of the 0^ state vanishes in the plot with reordering, and 
the point moves smoothly and approaches the exact en- 
ergy on the y-axis as usual [17]. These points are fitted 
by a first-order polynomial, which is shown as a dashed 
line. Unlike the fit without reordering, these points are 
on the line in a large range of the energy variance which 
can be used for the fitting. This makes the extrapolation 
procedure stable. 
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FIG. 3: (Color online) Energy vs. energy variance plot for the 
ground state and state of '^^Ge in /sp^g-shell. The open 
triangles are obtained by the CG method without reordering 
and the solid circles are with reordering. The inset shows a 
magnified view for the 
text for detail. 



-181 



energy around (A_ff ) ~ 0. See 



> 



>;-182 

0) 

c 
LU 



-183 







—m — 1 — ' — r- 


(a) / 

/ X ♦ * 


(b) 








A 






A 












A 


^^^^^^^^^^^^^^^ 




A 

A 

A 






A 

^ A 


exact j=6^i 













10 

Energy Variance (MeV^) 



20 20 40 
Number of basis states 



FIG. 4: (Color online) (a) Energy vs. energy variance plot 
of the ground state energy of ^^Ge in the /spgg-shell with 
the SCG wave function and reordering. Nm — 10 (orange 
crosses), 20 (green diamonds), 30 (red triangles), 40 (blue 
squares), and 50 (open black circles), (b) Extrapolated value 
vs. the number of basis states, Nm- The exact energy is also 
shown by the dotted line. The solid points and open blue tri- 
angles denote the extrapolated energies with first-order fitting 
and reordering, and those with second-order fitting without 
reordering, respectively. 



Figure |31Ja) shows the EV plot provided by the SCG 
method with a various number of basis states, Nm = 
10,20,30,40 and 50. At each N,n, the last 10 points 
are used to make the fitted line. While the extrapolated 
values of Nm = 10 apparently underestimate the ground 
state energy, the value converges as a function of Nm for 
Nm > 30. 

Discussing the stability of the extrapolated value is dif- 



ficult because all the energy- variance points with reorder- 
ing change by increasing the N. To discuss this stability, 
we show the extrapolated energy itself with the reorder- 
ing technique as a function of Nm- Figure S] (b) shows 
the convergence property of the extrapolated energy as a 
function of Nm- The extrapolated energy with reorder- 
ing using the last 10 points, or {Nm — 9)-th, (Nm — 8)-th, 
... iVm-th points, converges well within a few keV for 
Nm > 40. On the other hand, the extrapolated value 
without reordering shows a large fluctuation due to the 
kink in the EV plot. 

The inset of Fig. [3] shows the EV plot concerning 
J = O2 - The open triangles and the fitted line show 
the energy and variance without reordering. Its extrapo- 
lation underestimates 20 keV, the exact value, while the 
extrapolated value with the reordering technique agrees 
with the exact one within the 6 keV error. The reorder- 
ing technique reduces the error of the extrapolation even 
in the case of the simple EVE without reordering, and 
works well. 



C. ^*Ge in p/g9-shell 

In Ref. [13], we show a demonstration of the va- 
lidity of the MCSM and energy- variance by ^^Ge in 
the pf + (7g/2-shell model space, which consists of 
0/7/2,0/5/2,1^3/2,1^1/2 and O59/2 single-pM-ticle orbits. 
We use the PFG9B3 effective interaction [23| , which was 
also used in Refs. [13, HI]- The m-scheme dimension of 
the system reaches 1.7 x 10"'^'', which cannot be reached 
by the conventional Lanczos diagonalization technique. 
In the present work, this nucleus is taken as an example 
again. We show the results using the SCG method and 
by extrapolation with the reordering technique in Fig. [S] 
As you can see, the extrapolated energies with reorder- 
ing technique well converge in a relatively small number 
of basis states. In this case, the energies with reordering 
agree well with those without reordering, and those of 
the MCSM method which are shown in Ref.[l7|- 



V. SUMMARY 

We have discussed variational procedures, the SCG 
method, for nuclear shell-model calculations. Based on 
the SCG wave functions, we propose a new procedure of 
energy variance extrapolation, which can solve a complex 
problem due to shape coexistence. 

For example, the simple EVE is difficult to solve at the 
0^ state of ''^Gc with the JUN45 interaction. In this case, 
the anomalous kink appears in the EV plot, which shows 
the transition of the approximated wave function from 
the second lowest state to the lowest state. We discussed 
its origin in view of shape coexistence, and found that 
the reordering technique of the basis states allows us to 
extrapolate the eigenenergy successfully even in this case. 
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FIG. 5: (Color online) (a) Energy vs. energy variance plot 
for the O^"" and 2\ states of '^'*Ge in the p/gO-shell obtained 
by the SCG method, (b) Extrapolated energies with reorder- 
ing (filled circles) and without reordering (open triangles) as 
functions of A''™. 



We demonstrated that the reordering technique in the 
EV plot allows us to make a linear fit; and therefore, 



the reordering of the basis states makes the extrapola- 
tion procedure stable and suppresses the uncertainties 
of the extrapolated value. This procedure is expected 
to be quite useful into performing the precise estimate 
of nuclear energies based on large-scale shell-model cal- 
culations and no-core shell- model calculations [21|. For 
the estimatation of observables other than energy, such 
as quadrupole moment, the order which was used in the 
present work is not suitable for the extrapolation proce- 
dure. We are investigating a way to determine an order 
suitable for the extrapolation of these observables. 
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